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ABSTRACT 



Multifluid, Magnetohydrodynamic Shock Waves with Grain 
Dynamics II. Dust and the Critical Speed for C Shocks 

- - - Glenn E. Ciolek/'^ Wayne G. Roberge/'^ and Telemachos Ch. Mouschovias^ 

o 

O' cioleg@rpi.edu, roberw@rpi.edu 

Ph. 
< 

This is the second in a series of papers on the effects of dust on multifluid, 
^ ■ magnetohydrodynamic shock waves in weakly-ionized molecular gas. We investi- 

^ ■ gate the influence of dust on the critical shock speed, Vcrit-, above which C shocks 

C\l . cease to exist. Chernoff showed that fcrit cannot exceed the grain magnetosound 

^ ■ speed, V^ms, if dust grains are dynamically well coupled to the magnetic field. 

'nI" . Since Vg^s — 5 km s~^ in a dense cloud or core, the potential implications for 

models of shock emission are profound. We present numerical simulations of 
P_|. steady shocks where the grains may be well- or poorly coupled to the field. We 

Q . use a time-dependent, multifiuid MHD code that models the plasma as a system 

of interacting fiuids: neutral particles, ions, electrons, and various "dust fiuids" 



c/3 . 

. comprised of grains with different sizes and charges. Our simulations include 

grain inertia and grain charge fluctuations but to highlight the essential physics 
^ , we assume adiabatic flow, single-size grains, and neglect the effects of chemistry. 

We show that the existence of a phase speed does not necessarily mean that C 
shocks will form for all shock speeds Vg less than f^. When the grains are weakly 
coupled to the fleld, steady, adiabatic shocks resemble shocks with no dust: the 
transition to J type flow occurs at fcrit ~ 2.76V^a, where V^a is the neutral 
Alfven speed, and steady shocks with Vs > 2.76KiA are J shocks with magnetic 
precursors in the ion-electron fluid. When the grains are strongly coupled to the 
fleld, fcrit = min (2.76l^A, Vgms)- Shocks with Wcrit < "^s < ^ms have magnetic 
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precursors in the ion-electron-dust fluid. Shocks with Vg > Vgms have no magnetic 
precursor in any fluid. We present time-dependent calculations to study the for- 
mation of steady multifluid shocks. The dynamics differ qualitatively depending 
on whether or not the grains and field are well coupled. 

Subject headings: diffusion — dust, extinction — ISM: clouds — ISM: magnetic 
fields — MHD — plasmas — shock waves — waves 

1. Introduction 

Shock waves in weakly-ionized plasmas, such as the molecular outflows around young 
stellar objects, may have a multifluid structure, in which the charged and neutral compo- 
nents of the plasma behave as distinct fluids. If the shock speed is less than the speed with 
which the charged fluid communicates compressive disturbances, the charged particles are 
accelerated ahead of the neutral gas in a "magnetic precursor" (MuUan 1971; Draine 1980). 
The streaming of charged particles through the neutrals also accelerates, compresses and 
heats the neutral fluid. If it remains sufficiently cool, the neutral fluid is everywhere super- 
sonic in a frame at rest in the post-shock gas. In this case all hydrodynamical variables are 
continuous and the entire structure is termed a "C shock" (Draine 1980). Alternatively the 
neutral fluid must undergo a transonic transition. This may occur smoothly via continous 
flow through a sonic point (a "C* shock", see Chernoff 1987 and Roberge & Draine 1990, 
but note that Chieze, Pineau des Forets, & Flower 1998 question the reahty of C* shocks) 
or abruptly at a viscous jump front (a "J shock"). 

Simulations of the different types of multifluid MHD shocks reveal distinct observational 
signatures, including the strengths and velocity profiles of atomic and molecular emission 
hues (Chernoff, HoUenbach, & McKee 1982; Draine & Roberge 1982; Draine, Roberge, & 
Dalgarno 1983, henceforth DRD83; Smith & Brand 1990a,b; Smith 1991a,b; Smith, Brand, & 
Moorhouse 1991a,b; Kaufman & Neufeld 1996a,b; Timmermann 1996; Neufeld & Stone 1997; 
Pineau des Forets et al. 2001), the ortho/para ratio of H2 (Timmermann 1998, Wilgenbus 
et al. 2000), and the abundances of various chemical species in the post-shock gas (e.g., 
DRD83; Flower, Pineau des Forets, & Hartquist 1985; Pineau des Forets, Flower, & Dalgarno 
1988; Flower & Pineau des Forets 1994; Flower et al. 1996; Bcrgin, Neufeld, & Mclnick 
1998). When the simulations are used to model the rich array of spectral observations now 
available for molecular outflows, a combination of C- and J shocks is broadly satisfactory, but 
an acceptable explanation for important aspects of the observations remains elusive (e.g.. 
Smith, Eisloffel, & Davis 1998; Rosenthal, Bertoldi, & Drapatz 2000). Consequently, the 
theory of multifluid shocks remains of considerable interest. 
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In this paper we reexamine a fundamental prediction of the theory, namely, the "critical 
speed," Vcrit, above which C type solutions fail to exist. According to most prior work, 
the transition from C- to J-type solutions occurs when the neutral flow becomes too hot 
to remain subsonic. For shocks in dense molecular gas, the transition is caused either by 
the onset of H2 dissociation (which decreases radiative cooling of the neutral gas) or by 
self- ionization (which increases heating by ion-neutral scattering). The precise value of Wcrit 
is sensitive to the rates of various atomic and molecular processes. DRD83 estimated that 
49 > 

fcrit ^ 40 km s ^ for preshock densities in the range 10^ < < 10^ cm ^. Analogous 
calculations by Smith (1994) give similar values, though recent simulations by Le Bourlot et 
al. (2002) predict somewhat larger critical speeds (e.g., fcrit ~ 70 km s~^ for riu — lO^cm"^). 
In any case, shocks in dense molecular gas would be C type for speeds >40-50km s~^ if the 
critical speed was determined solely by the temperature of the neutrals. 

However Vcvit is also constrained by the "signal speed" for compressive waves in the 
charged fluid, since shocks faster than the signal speed must necessarily be J type. In a cold, 
weakly-ionized molecular gas, Vgig ~ B/^Anpc, where B is the magnetic field strength and 
Pc is the mass density of charged particles. If pc is taken to be the density of atomic and 
molecular ions, then Vsig = Via, where ViA ^ 1000 km s~^ is the ion Alfven speed. In this view 
the constraint imposed by the signal speed would be of no practical significance. However 
dust grains are also (mostly) charged (e.g., see Fig. 1 of Paper 1) and, to the extent that 
the grains are "loaded" onto the field hues, their contribution to pc also should be included 
(McKee, Chernoff, & HoUenbach 1984). For example, if pc is taken to be the total dust 
density, then Vgig — Vg^is, where V^ms ~ 5 km s~^ (see eq. [28]) is the grain magnetosound 
speed. A critical speed this small would have profound implications for models of shocked 
molecular outflows. Flower & Pineau des Forets (2003) have computed the abundance of 
charged grains in clouds including polycyclic aromatic hydrocarbons (PAHs) and cosmic-ray- 
induced UV photons. They suggest that only a small fraction of the total grain mass in a 
plasma containing PAHs is charged and coupled to the field, generally yielding magnetosound 
speeds ~ 50 km s~^. However, in the environments we wish to study, namely, dense, dark 
clouds and cores, it is likely that most PAHs have condensed or frozen onto the icy mantles of 
larger grains (Ehrenfreund & Charnley 2000; Bernstein et al. 2003; Gudipati & AUamondola 
2003; Kriigel 2003), and therefore should not affect the charge state of grains to such a large 
degree. 

The critical speed t'crit would equal the grain magnetosound speed T^ms only if all the 
grains are well coupled to the magnetic field fines (Chernoff 1987). However, the grains 
may be only weakly coupled (e.g., Pilipp et al. 1990; Draine & McKee 1993) and grain 
charge fluctuations limit the effective mass of the loaded grains to some fraction of the 
total dust mass (Ciolek & Mouschovias 1993). Consequently, detailed numerical simulations 
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are generally required to determine T;crit- Although dust has been included in many detailed 
studies of shocks, most have neglected the inertia of the grains (i.e., the momentum equation 
for the dust has been omitted, on the assumption that the dust motion is force free). Because 
propagating waves in the grain fluid are suppressed when the grain inertia is neglected, 
studies of this type cannot be used to examine the connection between grain loading and 

In this paper we present numerical calculations which include the grain inertia, and 
extend Chernoff 's work by also including partial grain-field coupling and grain charge fluctu- 
ations. Our plan is as follows: In § 2 wc summarize our assumptions and the free parameters 
in our models. The signal speed in a dusty plasma is discussed in § 3, where we examine 
the effects of dust on the phase velocities and damping rates of compressive waves. In § 4 
we present numerical solutions for steady and time-dependent shock waves. Our results are 
summarized in § 5. 

2. Modeling Assumptions and Parameters 

We simulate time-dependent, MHD shock waves propagating in weakly-ionized plasmas 
using a finite difference code which is described in detail elsewhere (Ciolek & Roberge 2001, 
henceforth Paper I). We consider plane-parallel shocks propagating along the ±z direction 

of a Cartesian frame {x, y, z). Fluid variables are assumed to depend only on z and time 
t. We restrict our discussion to "perpendicular shocks" with a magnetic fleld of the form 
B = B{z,t)x. Fluid motions in the x-direction (i.e., along the magnetic fleld) arc uncoupled 
from those in the y and 2;-directions and are ignored. Self gravity of the gas is also neglected. 

Our code models the plasma as a system of up to nine interacting fluids: neutral particles 
(denoted by subscript n), electrons (e), ions (i), large dust grains (g~, g", and g+) and very 
small grains or PAHs (sg^, sg°, sg+) with charges —1, 0, and +1 times the proton charge, e. 
The governing equations presented in Paper I arc solved here for the special case of vanishing 
PAH or very small grain abundances since, as noted in the preceding section, they are likely 
to be depleted onto large grains in dense cores. 

Our code has the capacity to include radiative cooling of the neutral gas by H2, CO, and 
H2O molecules, as well as cooling by collision-induced dissociation of H2. For the calculations 
reported in this paper, however, we "turned off" the neutral cooling. The assumption of 
adiabatic flow was adopted for two reasons: First, the structure of a steady adiabatic shock 
with constant fractional ionization can be calculated exactly (Chernoff 1987; Smith & Brand 
1990c), and this provides a valuable benchmark on our numerical results. Second, Chernoff 
(1987) showed that shocks in a weakly- ionized plasma undergo a C- to J-type transition 
at I'crit = 2.76\4a, where VnA — B/^/^np^ is the neutral Alfven speed. This conclusion 
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holds if (i) the charged particles are well coupled to the magnetic field and (ii) the shock 
is subAlfvenic in the charged fluid. Since these criteria are always satisfied by the ions, 
the occurrence of a C-J transition at any other critical speed, if found, can be attributed 
unambiguously to the effects of dust. With neutral cooling turned off, only the abundances 
of the major species H, H2, and He (on which the "source terms" for energy and momentum 
transfer between the charged and neutral fluids depend; see Paper I) are required. We omit 
atomic hydrogen since the shocks considered here are too slow to be dissociative and fix the 
other abundances at n(H2) = O.Stt-h and n(He) = O-lnn- 

The ion fluid is composed of generic atomic (subscript a"'') and molecular (m"'') ions 
with identical masses, ma+ = mra+ = 25amu. The ion abundances are described by a 
rudimentary chemical network which allows for changes in the fractional ionization caused 
by dissociative recombination, etc (see Paper I). In this paper we "turned off" the chemistry 
so that steady solutions could be calculated exactly. We neglected the inertia of the ions 
and electrons because they attain force-free motion on length scales that arc many orders of 
magnitude smaller than the flow scale, and are therefore impractical to track with our finite 
difference scheme. Thermal pressure forces on the ions and electrons are negligible compared 
to electromagnetic forces and gas drag, and so were ignored. Because the heat capacities of 
the ion and electron fluids are small, their temperatures were calculated by requiring these 
fluids to be in thermal balance (see Chernoff 1987). 

The large grains (henceforth simply "grains") were assumed to be spherical with iden- 
tical radii, Og, and to be composed of material with density 3gcm~^. Realistic simulations 
with a spectrum of grain sizes will be presented elsewhere. Our purpose here is to elucidate 
the essential physics. 

The free parameters in our calculations arc listed in Table 1. Parameters that describe 
the undisturbed preshock plasma have realistic values for a typical cloud core, except that 
the cosmic ray ionization rate has been artiflcially adjusted so that all calculations have 
the same fractional ionization. We considered two simple dust models, in which the dust is 
either weakly coupled (the WC model) or strongly coupled (the SC model) to the magnetic 
field. Both models have a reahstic dust-to-gas mass ratio, but the grain radii were chosen 
artificially to give large and small values for the grain Hall parameter. 



where Qg is the charged-grain gyrofrequency and r^f^g grain-neutral collision time. 

We give values of Fg and other pertinent quantities in Table 2. 
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3. Compressive Disturbances in a Weeikly Ionized, Dusty Plasma 

3.1. Wave Modes 

A weakly-ionized, dusty plasma supports a variety of waves with phase velocities and 
damping times that can be obtained by solving the appropriate dispersion relations. Pilipp et 
al. (1987, hereafter PHHM87) solved the dispersion relation for Alfven (i.e., noncompressive) 
waves including the effects of partial field coupling and charge fluctuations. PHHM87 found 
that, when the grains are well coupled to the field lines, their added inertia reduces the phase 
velocity relative to a dust-free plasma that is otherwise identical. Ciolek & Mouschovias 
(1989, hereafter CM89) carried out the analogous calculations for waves propagating at 
an arbitrary direction relative to the field. Here we consider the special case of CM89's 
analysis for the magnetosonic (i.e., compressive) modes relevant to the formation of magnetic 
precursors. Since some of the physics is both interesting and not published elsewhere, we 
describe the results in some detail. 

We obtained the dispersion relation by linearizing the equations of motion about a 
uniform, static unperturbed state and searching for monochromatic plane-wave solutions of 
the form 



for the first-order perturbation qi in variable q. The perturbations were assumed to have 
the geometry appropriate for magnetosonic modes. We linearized the equations of motion 
from Paper I with several modifications: We omitted the equations of motion for very small 
grains on the assumption they are depleted (see §2). We omitted the continuity equations 
for the ions and grains because, for the very low fractional ionizations of interest here, their 
density perturbations have virtually no affect on the dynamics. We included the momentum 
equation for the ions (the ion inertia is neglected in Paper I and in our numerical code) 
since otherwise waves in the ion fluid would be suppressed. We assumed for simplicity that 
the perturbations are adiabatic, and therefore omitted the energy equation for the neutral 
fluid. Finally, in writing the magnetic induction equation we assumed that magnetic flux 
is frozen into the electron fluid. This is more general than the corresponding assumption 
in Paper I (flux freezing in the ion fluid) though equivalent in the limit of vanishing small 
grain abundance. Our results for the WC and SC models are shown in Figures 1 and 2, 
respectively, where = Re[a;]/A; is the phase velocity and Tdamp = — 1/Im[a;] is the damping 
time. Several wave modes are apparent. 




(2) 
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3.1.1. Ion Magnetosound Waves and Quasiparticle Oscillations 

At very short wavelengths only the ions and electrons are coupled to the magnetic field. 
In the ion magnetosound mode (labeled "ims" in Fig. 1-2), the ions, electrons and magnetic 
field oscillate against a background of fixed neutrals and grains. The phase velocity is the 
ion magnetosound speed, 

Mms ^ Ma = (3) 

where l^ms — Via because we have neglected thermal pressure of the ions. For the core 
conditions adopted here, \4ms ~ 1000 km s~^ (Table 3). 

Ion magnetosound waves do not exist for wavelengths larger than some value Am^^. In 
the WC model, AmaS"* is the familiar cutoff imposed by ion-neutral friction, 

>^ = ^7iV,^.T^;^l (WC model) (4) 

(e.g., Kulsrud & Pearce 1969), where r^^^^ is the ion-neutral drag time.^ In the SC model, 
the cutoff for ion magnetosound waves is the electrostatic wavelength, 

A^r ) = Aeiect ^ (SC model) . (5) 

^ 'qp 

For A > Acicct the ion magnetosound mode transforms into a nonpropagating^ "quasiparticle" 
mode (labeled "qp" in Fig. 1-2). In the presence of charged grains, the ions and electrons 
collectively form a single "quasiparticle fluid" of electron-shielded ions with effective charge 



The ions and electrons gyrate about the field fines at a constant frequency 

= (!!iz^) fi, = (iv^) ft ^ fiv \ n,. (7) 



^In general, t^^^^ denotes the drag time for particles of type i to be slowed by elastic scattering with 
particles of type j, and T^j'^'^g*^ denotes the drag time for particles of type i to be slowed by scattering with 
all other particle types. 

^Note that d{logv^)/d{logX) = 1 on the branch labeled "qp". 

^The quasiparticle mode also appears in the WC model (Fig. 1). However, quasiparticle effects are 
insignificant in the WC model: the abundance of charged grains, and hence QeS, are « 60 times smaller 
than in the SC model (see Table 2). As a result, Oqp is so small that the quasiparticle oscillations are 
evanescent. 
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where Qi = eB/rriiC is the ion gyrofrequency. For detailed discussions on the physics of 
quasiparticle oscillations, see Ciolek & Mouschovias (1993, 1994, 1995). 

In both models, propagating ion magnetosound waves are damped by ion-neutral drag 
on a time scale 

(ims) _ r, (i,n) (s ^{ims)^ 

where the factor of 2 occurs because of equipartition between the kinetic and magnetic field 
energies. ForA> A^r^ ion- magnetosound waves are evanescent: the ions and magnetic field 
diffuse through the rest of the plasma with damping on the diffusion time scale, 

\2 

(ims) _ ^ (\ ^ \(ims)x /Qv 

^damp - (ims) > ^"^^ ^ 

TJ" I'ims^max 



3.1.2. Grain Magnetosound Waves 

At somewhat larger wavelengths, grains participate in the wave motions provided they 
are coupled to the magnetic field. Thus, grain magnetosound waves (labeled "gms" ) appear 
in the SC model but not in the WC model. In the former, g~ and g"*" grains are directly 
attached to the field lines and the neutral grains are effectively attached by charge fluctu- 
ations.^ Accordingly, ions, electrons and dust all contribute to the inertia of the wave and 
the phase velocity is the grain magnetosound speed, 

^gms = . . , f ^ y (10) 
V47r (pi + Pe + Pg) 

where pg is the total density of charged and uncharged dust. Because the dust density is 
much larger than the ion density, the grain magnetosound speed is much smaller than the 
ion magnetosound speed (V^gms ^ms, see Table 3). 

Grain magnetosound waves are allowed for wavelengths larger than 

^fems) ^ TrVg-^, 

^ 'qp ' drag ^ 'g 

where Vg-ms (— V^ms, since rig- ^ UgO, ng+ — see Table 2) is the magnetosound speed in 
the fluid of charged grains and Qg is the grain gyrofrequency. The minimum wavelength 



^Attachment of the neutral grains occurs if the time for a neutral grain to capture an electron or ion 
is much less than the g°-neutral drag time, T^£g"' (PHHM87; CM89; Ciolek & Mouschovias 1993). This 
inequality is satisfied by the SC and "intermediately-coupled" IC (see § 4.2.1) models. 
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A^™^^ signifies some interesting and unique plasma physics: propagation of magnetosound 
waves in the grain fluids is aided by an electrostatic attraction between the positively-charged 
quasiparticles and the negatively charged grains (CM89). For A > X^l^"\ the electrostatic 
attraction between the quasiparticles and charged grains allows the charged grains to be 
"dragged along" with the magnetically attached quasiparticles and the field (see Ciolek & 
Mouschovias 1993, § 3.1.2). The maximum wavelength for grain magnetosound waves is 
determined by gas-grain friction. The cutoff occurs at 



A 



(gms) 



4 T/ (g.n) 



1 + 



drag I 



(12) 



(CM89). For A > A^^Tx^ , the waves are evanescent and the combined plasma of ions, electrons 
and dust diffuses relative to neutrals. 

At small wavelengths (i.e., A ~ A^™^^), grain magnetosound waves are damped primarily 
by the drift of quasiparticles and magnetic field relative to the charged grains. The damping 
time is 

1 -(gms) ^^3^ 



.(gms) 



damp \(gms) y 



if A ~ A 



g~ms 



and this "damping branch" therefore has a slope rflogTdamp/c^log A = 2 in Fig. 26. At longer 
wavelengths the diffusion of quasiparticles with respect to grains is significantly lessened and 
the waves are damped by the frictional drag on the grains and ions exerted by the neutrals 
(which are essentially motionless at these frequencies) . In this regime the damping timescale 
becomes independent of A and the damping curve in Fig. 26 flattens. The damping rate is 
the sum of the damping rates for ion-neutral and grain-neutral friction, so that 



J. (gms) 
damp 



Pi/Pg 
(i,n) 



T, 



drag 



(g.n) 
drag 



if A 



(gms) 



< A < A 



(gms) 



where 



.(g,n) - X] 



Pg* 



(14) 



(15) 



'drag <5 Pg ''"drag 

and the factor of 2 again appears because of equipartition between the kinetic and magnetic 
field energies in the wave. For A > AmSf the (evanescent) waves are damped on the diffusion 
time, with 

ffA>ALT\ (16) 



.(gms) 



damp 



T/ \ (gms) 
•■^gms '^max 
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3.1.3. Neutral Sound and Magnetosound Waves 

Oscillations of the neutral fluid occur for all wavelengths. At small wavelengths, adi- 
abatic neutral sound waves (labeled "ns") occur wherein the neutrals oscillate in a sea of 
fixed ions and grains. The phase velocity is the adiabatic sound speed, 

c^ = y.c„.(^)"\ (17) 

where 7 is the adiabatic index of the gas (Table 3). Sound waves in the neutral fluid are 
damped by collisions between neutrals and the other fluids, so that 

-1 -1 

damp 

At longer wavelengths the ions and grains excite neutral magnetosound waves (labeled 
"nms") consisting of longitudinal compressions of the bulk fluid (i.e., the neutrals plus the 
electron-ion-grain plasma) and magnetic field. The phase velocity is the neutral magne- 
tosound speed, 

Vnms-{Cl + V^\y^\ (19) 

(Table 3). Since the neutral fluid density far exceeds the combined densities of all other fluids, 
the inertial response of the wave is almost entirely due to the neutrals. Thermal-pressure and 
magnetic forces — the latter affecting the neutral fluid indirectly through collisions between 
neutrals and charged particles — provide the restoring force required for wave propagation. 

Propagation of neutral magnetosound waves becomes possible for wavelengths greater 
than the minimum neutral magnetosound wavelength, 

ir I 'drag 
'nms ' 



aLT^ = ^ Ut^ r^S^ (20) 



(CM89). Note that for T4a ^ C'ad, (^nA/^nms) — Ka- The minimum neutral magnetosound 
wavelength is the perpendicular (k J- B) analog to the lower cutoff wavelength for neutral 
Alfven waves {k \\ B) in ion-neutral (Kulsrud & Pearce 1969) and dusty (PHHM87) plasmas. 
The damping of neutral magnetosound waves is due to ambipolar diffusion of the plasma 
and magnetic field with respect to the neutrals, with 

damp r,_2T/2 ^(n.tot) V ) 

'^nAMrag 

(see Fig. 16,26). For A > A^^^^ there also exists a non-propagating i\v^\ — 0) thermal- 
pressure-driven diffusion mode, in which the neutrals diffuse relative to a fixed background 
of plasma and magnetic field, which has Tdamp = A^/47r^Cn^T^rag°*^" 
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3.2. Propagation of Small- Amplitude Disturbances 

The preceding section shows that disturbances in a dusty plasma can excite a variety of 
waves, with phase velocities and damping rates that span several orders of magnitude. We 
are ultimately interested in how compressive waves in the charged fluids affect the formation 
of C shocks, and we do this in §4 by simulating large-amplitude disturbances (i.e., shocks). 
However, it is instructive to consider first a simpler problem, the propagation of small- 
amplitude, Gaussian wave packets, where interpretation of the dynamics in terms of different 
modes is much simpler: the narrow range of length scales in such a packet limits the number 
of modes that are excited, and the linearity of the disturbance allows us to interpret the 
results using §3.1. 

At time t — we superimposed a perturbation 

{z - Z/2f 



B{z) = 5B exp 



02 



(22) 



on a static, uniform initial state corresponding to either the WC or SC model. We centered 
the pulse initially on a computational domain [0, Z] and followed its evolution with our 
numerical code. The code uses transmissive boundary conditions, which allow the pulse to 
exit the domain freely, i.e., without generating reflected pulses at the boundaries. Thus, the 
pulse should evolve exactly like an identical pulse in an infinite computational domain. In 
all cases we set 5S /So = 0.1 but chose a few different values of the packet width, Ib- 

Simulations of a wave packet with Ib — 4x10^"^ cm are described in Figure 3 for the WC 
model, where the left-hand panel shows fluid velocities and the right-hand panel shows the 
density and magnetic field. The upper and lower boxes show variables at times ti — 8.5 yr 
and t2 — 17 yr, respectively. These times are so short {ti,t2 <C T^rag, see Table 2) that 
the neutrals are virtually stationary. The grains, which receive their "marching orders" in 
the WC model mainly by colliding with the neutrals, are also motionless. The ions are not 
stationary, but the initial disturbance does not propagate: the "wave" packet simply sits at 
the center of the grid and decays away. The ions in Fig. 3 are diffusing through the grains 
and neutrals at the terminal drift speed determined by the balance between ion-neutral 
drag and the driving magnetic pressure gradient. The absence of wavelike (i.e., oscillatory) 
motions may seem surprising, since propagating ion magnetosound waves are permitted by 
the dispersion relation on the length scales ~ in the initial pulse (Fig. la). However the 
damping time for propagating waves in the ions is so short (~ 0.01 yr, see Fig. 16) that such 
waves have virtually no effect on the dynamics.^ We conclude that, on time scales of practical 



^Of course, ion magnetosound waves are suppressed by our numerical code, where the ion inertia is 
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interest, ion motions in a compressive disturbance are the diffusive motions associated with 
evanescent waves. 

Figure 4 describes the same initial disturbance as Figure 3 but for the SC model where 
grain magnetosound waves are allowed for wavelengths ~ £b (Fig. 2a). In contrast to the 
preceding example, the initial disturbance gives rise to two symmetric wave packets propa- 
gating in the ±z directions. That the pulse in Fig. 4 is indeed a packet of grain magnetosound 
waves can be confirmed by computing the propagation speed: the centers of the forward- 
and backward-propagating pulses travel a distance ~ 1.2 x lO"*^^ cm between ti = 8.5 yr 
(Fig. 46) and t2 = 77 yr (Fig. Ad). The propagation speed is therefore 5.0 km s~^, in ex- 
cellent agreement with V^^s (Table 3). The packet is also damped at the expected rate. 
Taking the width of the pulses (~ 2.6 x 10^^ cm) at time ^2 as a characteristic wavelength, 
we find from Fig. 2b that the damping time should be ~ 92 yr. Since the initial amplitude of 
each counterpropagating pulse is O.OSSq, the value ol B/Bq at the peak of the pulse should 
therefore be ~ 1 -|- 0.05 exp(— 77/90) = 1.022 at t2- This is very close to the numerical value, 
1.024, for the data in Figure Ad. Finally, the linear analysis for grain magnetosound waves 
predicts that the magnitude of the grain velocity^ should obey the relation 



in the limit of weak damping, where Cn is the isothermal sound speed (Table 3). At points 
near the peaks of the pulses, the numerical data for Fig. 4c and Ad satisfy this relation to 
better than one percent. 

Figure 4 also illustrates the different responses of the various fluids to the compressive 
disturbance. For example, the neutrals are unable to respond at the frequencies correspond- 
ing to length scales ~ £b and act simply as a stationary background. The grain velocity 
lags somewhat behind the ion velocity (Fig. 4a, c) because the time for the charged grains to 
respond to the magnetic field variations (Jlj^=63 yr, see Table 2) is not much smaller than 
the flow time scale (< lO^^cm/5 km s~^ = 65 yr). In fact the grains repond somewhat faster 
than this because of their attraction to the electron-shielded ion quasiparticles, which are 
attached to the magnetic field. Thus, the ion and grain motions are becoming synchronous 
by time ^2 = 72 yr (Fig. Ac). 



neglected. The smallness of the damping time compared to other dynamical time scales of interest justifies 
this omission. 

^The neutral and charged grains behave as a single fluid in our simulations of wave packets because the 
charge fluctuation time is small compared to the flow time scale {'^ is/v). 




(23) 
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Figure 5 describes another Gaussian packet in the SC model, but with a larger width, 
Ib — 1.1 X 10^^, that is close to the maximum wavelength for propagating grain magnetosound 
waves (Fig. la). As expected, the motion is largely diffusion of the ion-grain plasma with 
respect to the neutrals, with some additional effects caused by the presence of Fourier com- 
ponents at wavelengths short enough for propagating grain magnetosound waves. Thus, the 
magnetic field pulse is more smeared out than the analogous pulse in Fig. 3, and has "rounded 
shoulders" (compare Figs. 3d and bd). The damping time scale (~ 330 yr, see Fig. 16) for 
the wave modes is much less than the decay time (> 5 x 10^ yr) for the longer-wavelength 
diffusion modes. Grain magnetosound waves are therefore overwhelmingly damped before 
wave packets can "break out" of the initial Gaussian pulse. After times comparable to the 
response time of the neutrals (~ T^ra|'', see Table 2) the neutrals acquire small nonzero veloc- 
ities (Fig. 5c) away from the center of the initial pulse, where a density depression appears 
(Fig. bd). 

Figure 6 describes an even wider packet with = 2 x 10^^ cm in the SC model. The 
width now exceeds the minimum wavelength for (propagating) neutral magnetosonic waves 
(Fig. la). At times shorter than ~ T^ral'' (Table 3), the velocity of the neutrals lags that of 
the ion-grain fluid (Fig. 6a, 6). At much later times (Fig. 6c,(i), the ions, grains, and neutrals 
behave as a single fluid and packets of neutral magnetosound waves start to propagate in 
the ±2; directions. Because the neutrals are effectively loaded onto the magnetic field hnes 
on these scales, the signal speed is greatly reduced. The data for Fig. 6 imply a propagation 
speed of about 0.57 km s"^, close to the neutral magnetosound speed (Table 3). 



4. Simulations of Multifiuid Shock Waves with Dust 

We now consider shock waves, which were simulated as follows. At an initial time, t = 0, 
the flow was assumed to be a J shock propagating in the -\-z direction, with the jump front 
located at the left-hand boundary of the domain, z — The preshock [z > 0) region was 
assumed to be uniform and static initially with values of the density, magnetic field strength 
and other variables taken from either the WC or SC model. Initial values of the post-shock 
{z < 0) variables were calculated by choosing some initial shock speed, fg i, and applying 
the adiabatic jump conditions (assuming a frozen-in magnetic field) aX z = 0. The fiow at 
later times was determined using our numerical code. Transmissive boundary conditions 
were applied at both boundaries to simulate a domain extending to z = ±00. Adiabatic fiow 
was assumed and the chemistry was turned off; the effects of cooling and chemistry will be 
discussed in later papers. 
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4.1. Dust and the Structure of Steady Shocks 

We computed a sequence of steady solutions for increasing shock speed, Vg, by integrating 
the equations of motion to a steady state. For each speed we computed solutions for both 
the WC and SC models in order to compare effects of dust in the two cases. Instead of 
the shock speed, it is useful to characterise the steady solutions by the neutral Alfven Mach 
number, 

M„A ^ (24) 

or the grain Alfven Mach number, 

MgA ^ (25) 

where V^a = Vgms is the grain Alfven speed. Recall that for a two-fluid (ion+neutral) plasma 
adiabatic shocks have Vcrit = 2.76V^a (Chernoff 1987) or Vcrit — 1.38 km s~^ for the Alfven 
speed adopted here (Table 3). To the extent that our solutions resemble the two-fluid case, 
we expect C type solutions only for shock speeds much smaller than the speeds (> 10 km s~^) 
in shocked molecular outflows. Thus, our numerical results bear little resemblance to real 
shocks: dust physics is the issue of interest here. 

In fact we expect the WC solutions to be virtually identical to two-fluid shocks with no 
dust: except for regions very close to a jump front, the grains should move with the neutrals, 
and therefore have no role in momentum or energy transfer between the charged and neutral 
fluids. This is verified by Figure 7. The symbols show velocities computed numerically for 
a shock with M^a = 2.1 using preshock conditions from the WC model. The smooth curves 
are the exact solution for a steady, two-fluid shock of the same speed (Chernoff 1987; Smith 
& Brand 1990c). 

Figures 8-11 describe a sequence of steady shocks with Alfven Mach numbers increasing 
from MnA = 2.1 to M^a = 5.0. It is difficult to pinpoint the C-J transition in these solutions 
because J-shocks are weak near the transition point (the compression ratio of the shock is 

unity precisely at the C-J transition) and because artificial viscosity smears out the J front. 
However it is apparent from Figs. 8-11 that (i) the transition is at or near MnA = 2.76 
for the WC model; and (ii) the transition occurs at or near the same M^a value for the 
SC model. The first result is expected for reasons noted above. The second result occurs 
because the critical speed, Wcrit — 1.38 km s~\ is smaller than the grain magnetosound speed, 
Vgms — 5 km s~^. That is, the transition to J- type flow occurs because the neutrals become 
too hot, not because the shock speed exceeds any signal speed. 

The principal difference between the WC and SC solutions in Fig. 8-11 is that the 
magnetic precursor is thinner in the SC models. Draine (1980) pointed out that the width of 
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the precursor is is inversely proportional to the rate per unit volume, F, at which momentum 
is exchanged between the charged and neutral fluids. In the WC models momentum transfer 
is due almost entirely to ion-neutral scattering. Dust contributes virtually nothing because 
the grains are almost at rest in the neutral fluid. However, in the SC model grain-neutral 
scattering actually contributes somewhat more to F than ion-neutral scattering (note that 
''"drag^ ^ '^drag Table 2), and this makes the precursors thinner. 

We conclude that the signal speed in the grain fluid has no effect on the C-J transition 
for adiabatic shocks if Vgms > 2.76MnA. Wc say that C shocks are "cooling limited" in this 
regime. To explore the regime Vgms < 2.76A'/„a! wc constructed a "heavily loaded" (HL) 
grain model identical to the SC model (i.e., the grains are well coupled to the field) but with 
an extremely large dust-to-gas mass ratio, Pg/pn = 0.5, chosen to make V^mg = 0.71km s"^. 
Figure 12 shows a sequence of steady shocks with increasing speeds computed for the HL 
model. The solutions are C type for MgA < 1 and J type for MgA > 1, as predicted by 
Chernoff (1985). We say that C shocks are "signal limited" in this regime, since the C-J 
transition occurs when Vs exceeds the signal speed (=V^ms) in the charged fluid. Notice that 
the ions are forced to move with the grains in these models, since both are coupled to the 
magnetic field. Consequently there is no magnetic precursor in any fiuid for shocks with 
MgA > 1' 



4.2. Dust and the Formation of Magnetic Precursors 

Our study of small-amplitude disturbances shows that the response of charged particles 
to a short- wavelength disturbance depends dramatically on whether or not the dust is coupled 
to the magnetic field. It is therefore of interest to study how a magnetic precursor forms 
in these two different cases. Figures 13 through 15 compare shocks in the WC and SC 
models for a sequence of increasing times. Except for the grain models, both calculations 
started with identical initial conditions: J shocks with identical compression ratios at z — 
and uniform, static plasma for 2; > 0. The steady solutions for the SC (Fig. 156) and WC 
(Fig. 15c) models are similar, consistent with the fact that the steady state has a shock speed 
Vs < 2.76Ka (see §4.1). However, the intermediate solutions are dramatically different. 

Our hnear analysis (§3.1) suggests that differences in the ffows might be caused by the 
different mechanisms for transmitting compressive disturbances in the two cases. Strictly 
speaking the disturbances in Fig. 13-15 are nonlinear; however the linear analysis should give 



^Thc slight offsets in Fig. 12c, d between the velocities of charged and neutral particles is a numerical 
artifact: they are only one grid zone wide. 
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a reasonable description of the charged fluids, since Vi <^ Vims and Vg <^ ^ms- We therefore 
adopt the view that the discontinuity imposed by the initial J shock "excites" modes with 
all possible wavelengths initially, and interpret the subsequent evolution of the flow in terms 
of these different modes. 

Consider first the WC solution. The earliest times shown (Fig. 13a, c and Fig. 14a) 
are much smaller than the dynamical response times of the dust (f^g^ ~ 4100 yr) and the 
neutrals {t^^^^^^ ~ 10"^ yr) so only the modes corresponding to pure ion motions are initially 
of interest. Since these times are also much longer than the damping time for propagating 
ion magnetosound waves (which are suppressed by our code in any case), the rapid ion 
motions in Fig. 13a, c correspond to evanescent ion magnetosound waves. To see that this 
interpretation is correct, recall (§3.1) that evanescent waves correspond to diffusive motions. 
A linear analysis of the ion equations of motion (Appendix A) shows that the diffusion 
coefficient for the ions is 



After a time t we therefore expect the "front" separating accelerated from stationary ions to 
have propagated a distance 



upstream from the initial J shock. In fact, the predictions of eq. (27) are in remarkably good 
agreement with the numerical solution for the WC model (Fig. 16). We conclude that a 
magnetic precursor forms via ion diffusion when dust is dynamically unimportant. 

For early times (Fig. 13a, c), the ions simply diffuse through an almost-stationary back- 
ground of neutrals and dust. After a few hundred years, these ion motions have com- 
pressed the magnetic field sufficiently to make the dynamical response time of the grains 
(~ ^ oc B~^) less than the flow time. At this point (Fig. 14a) betatron acceleration 
(Spitzer 1976) starts to drive the grain velocity toward the ion velocity. After another few 
thousand years the neutrals have responded significantly (Fig. 14c). Thereafter (Fig. 15a, c) 
the evolution slows as the charged and neutral fluids approach mechanical equilibrium, with 
magnetic pressure balanced by neutral drag in the charged fluid and ion-neutral drag bal- 
anced by thermal pressure in the neutral fluid. 

Our conclusion — that a magnetic precursor forms via ion diffusion rather than wave 
propagation when the dust is dynamically unimportant — has important implications. Ev- 
idently the conventional wisdom, that the ion magnetosound speed sets an upper limit on 
Vcrit when neutral cooling is strong, is incorrect. The upper limit in this situation must 
occur when the diffusion time scale becomes smaller than the time for ion-neutral scattering 
to accelerate the neutral gas. An important objective of our future work is to determine 




(26) 



-^-front 



(27) 
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whether this upper limit constrains real shocks with a spectrum of dust sizes and molecular 
coohng. 

Now consider the SC solution. Comparing 2;fi.ont(^) for the SC and WC models (Fig. 16), 
we see that the dynamics of the magnetic precursor are qualitatively different in the two 
cases. At very early times there is evidence for "i^/^ behavior" in the SC model (e.g., the 
curve in Fig. 16 is concave down for small times and the slope is very large near t = 0). We 
attribute this early time behavior to diffusive motions in the ion fluid, which is not strongly 
coupled to the grains for t ^ ^- At later times, -Zfront ~ Vg^st is a good approximation for 
the SC model. Clearly, the magnetic precursor forms in the SC model via the propagation of 
grain magnetosound waves. Now there is a well-defined signal speed, V^gms) and it is correct 
to say that Vcrit = Vgms- Oi course, this is precisely what we found in Fig. 12. 

It is of interest to note that a steady C shock forms at an earlier time in the SC model 
(pa 5 X 10^ yr. Fig. 15a) than in the WC model (pa 2.6 x 10^ yr. Fig. 15c). This can be 
understood for the following reasons: First, in the SC model, ions and grains move together 
and collide in parallel (i.e., simultaneously) with the neutrals. This gives a net collision or 
momentum-exchange time r^^^^c = ^d^f (1 + '^drag/i|V^ (e-g-, see eq. [18]). For the WC 
model, only coUisional drag between ions and neutrals is important, and T^rag°wc ~ "^drag- 
Inserting the values hsted in Table 2, we have 'rdr4°sW'^drag°wc ~ O-^- Hence, the coUisional 
effect of the ions and grains acting in concert results in momentum transfer that more rapidly 
"sweeps up" the pre-shocked neutrals in the SC model compared to that which occurs in the 
WC model. Secondly, the driving mean magnetic field gradient in the precursor region in the 
SC model is greater than in the WC model. As discussed above, this is due to the difference 
in the precursor formation mechanism between these models. Rapid ion diffusion disperses 
the magnetic field ahead of the shock in the WC model over a much larger region than in 
the SC model (Figs. 15a, c). Because the mean field gradient on the plasma is greater in the 
SC model than in the WC model, the neutrals thereby also experience a greater accelerating 
frictional force due to the larger gradient in plasma velocities occurring over a much shorter 
lengthscale. Thus, taken together, the combined effects of greater momentum transfer by 
collisions in a more narrowly confined magnetic precursor allows a steady C shock to form 
at a significantly earlier time in the SC model than in the WC model. 

4.2.1. Discussion 

For the SC model, which has a population of single-size grains with radius Og = 0.05 fim 
and Hall parameter Fg = 4.2, the grain magnetosound speed (eq. [10]) has the typical value 

1/2 

km s-^ (28) 



•'gms — 



B 



50 //G 



2 X 10 1 cur 



10 



-2 



.(Pg/Pn). 
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However, the single-size grain model that we have assumed here is an idealization to the 
interstellar medium; in actuality, grains have a distribution in grain sizes, normally spanning 
more than several orders of magnitude in radii. The very large grains in such a distribution 
will have Fg < 1 and will therefore not be coupled to the magnetic field (as in the WC 
model), while the smaller grains will instead have Fg > 1 and will be attached to the field 
(as in the SC model). The relevant question is then, what is the correct dust abundance 
that should be used to calculate the numerical value of T^ms when there is a distribution in 
grain sizes? 

A likely answer to this question is that one should use the fraction of the charged dust 
^ in the distribution that is coupled to the magnetic field. Neglecting grain-quasiparticle 
interactions, this corresponds to grains with Fg > 1. That the criterion Fg ~ 1 is the proper 
one for determining whether charged grains are attached to the magnetic field can be seen 
in Figures 17a, 6, which show C shock results for two "intermediately coupled" models, IC-1 
and IC-2, respectively. The two models have the same physical parameters as in the WC 
and SC models (including ion abundances), except that they contain different sized grains: 
IC-1 has grains with Og = 0.125 /im and Fg = 1.1, while IC-2 has grains with = 0.09 
and Fg = 1.3. The initial conditions for the shock flow are the same as for the SC and WC 
models in Figs. 13 - 15. As in the other models, a steady C shock eventually forms: in IC-1 
this occurs at 1.70 x 10^ yr, while in IC-2 it is at 1.25 x 10^ yr. Comparison of these models 
with the WC and SC models (Figs. 15a, c) shows that, even though the coupling of the 
grains to the magnetic field in the IC-1 and IC-2 models is not perfect, the shock structure 
has been significantly altered from that which occurs in the uncoupled WC model. Notably, 
the grains arc moving with the ions and field throughout most of the precursor region in the 
IC models, and the width of the precursor is substantially smaller than in the WC model. 
From this, we suggest that the criterion Fg = 1 can reasonably be used to separate grains 
with a continuous distribution of sizes into two distinct sub-populations of coupled (Fg > 1) 
and uncoupled (Fg < 1) charged grain fluids. 

The demarcation between these two sub-populations occurs at the maximum radius of 
flux-freezing in the grains, Og^fr, which is defined as the radius Og for which Fg = 1. From 



^It is assumed here that a non- negligible fraction of the grains are charged (see Table 2) ; this would also 
encompass grain models containing a population of very small grains for densities < lO'' cm~'^ (e.g., Nishi, 
Nakano, & Umebayashi 1991; see, also Fig. 1 of Paper I). As noted in § 1, rapid charge fluctuation and 
momentum transport between neutral and charged grain fluids can effectively couple even the neutral grains 
to the magnetic field (PHHM87; CM89; Ciolek & Mouschovias 1993). If, however, there is still an extant, 
non- negligible population of PAHs at these densities, the overall fraction of charged grains may be greatly 
reduced, with a corresponding increase in Vgms (Flower & Pineau des Forets (2003). 
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equation (1) it follows that 

a„fr = 1.0xlO-M cm. (29 

Grains with Ug < Cg^fr have Fg > 1; to calculate the mass fraction of grains contained in 
this sub-population, we use the well-known MRN grain distribution function, /(ag) oc flg"^'^ 
(Mathis, Rumpl, and Nordsieck 1977). Using this distribution function, the mass fraction of 
grains with Og < Og & is 



Pg{cig ^ Ctg,fr) / Q-gfi- 



1/2 



1 (flg,min/flg,fr) 



(^g,min/c^g,max) 



1/2 



(30) 



where ag^min and Cg^max are the minimum and maximum grain radii of the distribution. 
For our cloud/core models we may take Og^min = 0.01 fj,m, ag max = 0.3 /im, which yields 
Pg(ag < Og^fr) ~ 0.48pg. Hence, if we replace pg with Pg(ag < flgir) to estimate the grain 
magnetosound speed (eq. [10]) in interstellar plasma with a distribution of grain sizes, V^ms 
is only slightly greater than that predicted by equation (28), ~ 7 km s"^ — which is still 
considerably lower than the values of fcrit determined in the previously noted studies of 
"cooling limited" shocks. It is possible, then, that the existence and effects of a "signal- 
limiting" speed V^ms on the formation of shocks and magnetic precursors described in this 
paper will not only be qualitatively but also quantitatively correct even when there is a 
spectrum of grain radii. As stated above, we hope to definitively resolve this particular issue 
in future work. 



5. Summary 

The principal results of this paper are as follows: 

1. We calculated the dispersion relations for compressive waves in a dusty plasma for two 
models, where the dust is weakly coupled and strongly coupled to the magnetic field, 
respectively. In the "weakly coupled" (WC) model, the only propagating modes in 
the charged fluids are ion magnetosound waves at very short wavelengths and neutral 
magnetosound waves at very long wavelengths. In the "strongly coupled" (SC) model, 
grain magnetosound waves appear at intermediate wavelengths. 

2. We simulated the evolution of small-amplitude, Gaussian wave packets to see how dif- 
ferent plasma components respond to disturbances on various length scales. When the 
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length scales are small enough for decoupling of the charged and neutral particles to 
occur, the dynamics are dramatically different in the WC and SC models. In the WC 
model there is no wave propagation as such. Propagating ion magnetosound waves 
are allowed but strongly damped, and so have virtually no affect on the dynamics. 
Rapid diffusion of ions and electrons through the neutral gas and dust occurs, and the 
disturbance decays without propagating. In the SC model, packets of grain magne- 
tosound waves appear, and the disturbance propagates (with weak damping) at the 
grain magnetosound speed. Because the ions and electrons are tied to the magnetic 
field, the coupling of grains and magnetic field in the SC model suppresses diffusive 
ion motions. 

3. We simulated steady, adiabatic shock waves to compare the effects of dust in the WC 
and SC models. The WC model gives shocks that are virtually identical to shocks with 
no dust: steady shocks are C type for Alfven Mach numbers MnA < 2.76 and J type 
with magnetic precursors in the ion-electron fluid for MnA > 2.76. We say that steady, 
adiabatic C shocks are "cooling limited" in the WC model because the C-J transition 
occurs when the neutrals become too hot to remain everywhere supersonic. 

4. Steady, adiabatic C shocks are also coohng hmited in the SC model if Vg^ns > 2.76T4a- 
Their structures resemble their WC counterparts, with differences attributable to ad- 
ditional momentum transfer caused by grain-neutral scattering. 

5. In the SC model steady, adiabatic C shocks are "signal limited" if V^gms < 2.76T4a- 
In this case the C-J transition occurs when the shock speed exceeds the grain magne- 
tosound speed, as predicted by Chernoff (1985). Shocks faster than the grain magne- 
tosound speed have no magnetic precursor in any fluid. 

6. We simulated time-dependent shocks to see how magnetic precursors form in the WC 
and SC models. Consistent with our study of small-amplitude disturbances, the dy- 
namics are dramatically different in the two cases. In the WC case, a magnetic pre- 
cursor forms via diffusion of the ion-electron plasma through the dust and neutrals. In 
the SC case, the precursor forms via the propagation of grain magnetosound waves. 

Our calculations incorporated some very unreahstic assumptions in order to isolate the es- 
sential physics: Real grains have a spectrum of sizes and real shocks have chemistry and 
radiative cooling. Having understood the physics, our next investigation will determine 
realistic values of Vcnt by relaxing these assumptions. 

This work was supported by the New York Center for Studies on the Origins of Life 
(NSCORT) and the Department of Physics, Applied Physics, and Astronomy at Rensselaer 



Polytechnic Institute, under NASA grant NAG 5-7589. Comments from an anonymous 
referee are acknowledged with appreciation. 



A. Diffusion Coefficient for Ions in the Weakly Coupled (WC) Model 

The ion-magnetic field diffusion coefficient can be most easily obtained for this particular 
model by noting that the diffusion mode has the following characteristics (1) it has effective 
balance between magnetic forces and ion-neutral coUisional drag, and (2) the motions of 
the neutral and grain fluids are negligible compared to that of the ions. Moreover, as 
previously noted, we may also (3) neglect the effects of quasiparticlcs in the WC model. 
With these simplifying assumptions, the linearized versions of the plasma force (per unit 
volume) equation and magnetic induction equation are then, respectively, 

= -^f-fA^„, (Al) 

\ drag / 

m = -"'ai '*^> 

(see eqs. [17] and [9] of Paper I). Solving for Vi from equation (Al) and substituting that 
into (A2) yields the linearized magnetic field diffusion equation. 



3« ^ «^ ,<-.^ ^ . (A3) 



dt \47rpi '^'""Vo ' 

Fourier- analyzing this equation corresponds to the replacement 8/ dt — > —iu, and d/dz ik. 
Finally, using Tdamp = -1/Im[a;], A = 27i/k, Vi^s ^ Via = B^^o/ y/Airp^, and equation (4) 
yields the ion diffusion timescale (9). This can also be written as 

{' ) 

"^damp — ' (^4) 

^ion 

where Dion is the ion-magnetic field ambipolar diffusion coefficient (26). 
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Table 1. Model Parameters 



Parameter 


Symbol 


WC Model 


SC Model 


Preshock density (cm"^) 




2 X 10^ 


2 X 10^ 


Preshock magnetic field (/^G) 


Bo 


50 


50 


Preshock temperature (K)^ 


TnO 


10 


10 


Cosmic ray ionization rate'' (s^^) 


CcR 


8.4 X 10-19 


1.0 X 10-1^ 


Dust-to-gas mass ratio 


(Pg/Pn)o 


10-2 


10-2 


Grain radius (/xm) 




0.2 


0.05 



^Assigned arbitrarily, i.e., not determined by requiring the preshock gas to be in thermal 
balance. Assumed the same for all fluids. 

Values selected so that each model has the same degree of ionization (ni/nn)Q = 3.1 x 10~^. 
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Table 2. Properties of the Preshock Plasma 



Quantity 

Til J 


Symbol 


WC Model 


SC Model 


Ion-neutral drag time (yr) 


drag 


1.5 X 10^2 


1.5 X 10-2 


Neutral-ion drag time (yr) 


drag 


4.5 X 10^ 


4.5 X 10^ 


Neutral-grain drag time (yr) 


drag 


1.1 X 10^ 


3.0 X 10^ 


Grain-neutral drag time (yr) 


drag 


1.1 X 10=^ 


3.0 X 102 


Inverse grain gyrofrequency (yr) 




4083 


63 


Grain Hall parameter 


Tg 


0.26 


4.2 


g" abundance^ 


Xg- 


3.7 X 10-^3 


2.2 X 10-^^ 


g° abundance 




2.2 X 10-^^ 


2.7 X 10-^2 


g"*" abundance 




8.0 X 10"^^ 


4.8 X 10-15 


Inverse quasiparticle frequency (yr) 
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^Fractional abundance relative to neutral particles. 



-26- 



Table 3. Characteristic Speeds 



Quantity Symbol Value (km s ^) 

Ion magnetosound speed Vims 860 

Grain magnetosound speed l^ms 5.0 

Neutral magnetosound speed Km s 0.56 

Neutral Alfven speed Ka 0.50 

Neutral adiabatic sound speed Cad 0.25 

Neutral isothermal sound speed Cn 0.19 
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Fig. 1. — Dispersion relation for the WC dust model, where the grains are weakly coupled to 
the magnetic field, (a) Phase velocity = uj^/kys. wavelength, A, where uj-e = Re[uj] and k — 
27r/A. Various wave modes are labeled. Propagating modes include ion magnetosound (ims), 
neutral sound (ns), grain magnetosound (gms), and neutral magnetosound (nms) waves. The 
quasiparticle (qp) mode is a nonpropagating oscillation (sec text). (6) Damping time Tdamp = 
— l/uji, where cui = Im[a;]. Modes with difi^usive damping have a slope (i(logrdamp)/c^(log A) = 
2. 
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Fig. 2. — Similar to Fig. 1 but for the SC dust model, where the grains are strongly coupled 
to the magnetic field. 
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01234 01234 
z (10^^ cm) z (10^^ cm) 

Fig. 3. — Temporal evolution of a Gaussian wavepacket in the WC model. The initial packet 
width and amplitude are £b = 4 x 10^"^ cm and 5B = O.IBq, respectively. For visual clarity, 
only every fourth data point is plotted, (a) Velocities of the ions {filled diamonds), dust 
grains (crosses), and neutrals (open circles) at time ti — 8.5 yr. Velocities are in units of 
the isothermal sound speed of the undisturbed state (= 0.19 km s~^). (b) Magnetic field 
(filled diamonds) and neutral density pn (open circles) at time ti, normalized to values for 
the initial undisturbed state, (c) Same as (a) but at a later time t2 — 17 yr. (d) Same as 
(b), but at time ^2- 
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Fig. 4. — Same initial conditions as Fig. 3 but computed for the SC model, (a) Velocities 
at time ti = 8.5 yr. (6) Magnetic field and density at time ti. (c) Same as (a) but at time 
h = 77 yr. (d) Same as (6) but at time t2- 
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Fig. 5. — Same as Fig. 4 but for a larger initial packet width: £b — 1.1 x 10^^ cm. (a) and 
(b) are at ti — 340 yr, while (c) and (d) are at t2 — 4410 yr. 
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Fig. 6. — Same as Fig. 5 but for an even larger width: Ib — 2x 10^^ cm. (a) and (6) are at 
= 1.2 X 10^ yr, while (c) and {d) are at ^2 = 1-7 x 10^ yr. 
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Fig. 7. — Symbols: Numerical predictions of the ion velocity {filled diamonds) and neutral 
velocity {open circles) in a steady shock computed for the WC model. Smooth curves: exact 
solutions for a steady shock with the same speed. The neutral and grain Alfven Mach 
numbers are indicated. 
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Fig. 8. — Comparison of the velocity profiles for steady, adiabatic shocks with = 2.1 for 
the WC (Fig. 8a) and SC (Fig. 86) models. The grain Alfven Mach number is also indicated. 
Both solutions are C shocks. 




Fig. 9. — Same as Fig. 8 but for M-^a — 2.8. Both solutions are weak J shocks. The C-J 
transition has occured because the neutral fluid is too hot to remain subsonic. 
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Fig. 10. — Same as Fig. 9 but for M„a = 3.8. In both cases, the solutions are J shocks with 
magnetic precursors in the charged fluid. In the SC model, the grains move with the ions 
and electrons. In the WC model, the grains are nearly at rest in the neutral fluid, with a 
small velocity difference caused by betatron acceleration in the magnetic precursor. 
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Fig. 11.— Same as Fig. 8 but for M^a = 5.0. 
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Fig. 12. — Velocity profiles for steady shocks in the "HL" model, where the dust is well 

coupled to the magnetic field and the grain magnetosound speed is 0.72 km s~^. In this case 
the C-J transition occurs at MgA = 1, where the shock speed exceeds the signal speed in 
the charged fluid. Shocks with MgA < 1 (Fig. 12a, h) are C shocks with magnetic precursors 
in the charged fluid. Shocks with AfgA > 1 (Fig. 12c, d) are J shocks with no magnetic 
precursor in any fluid. (The slight velocity offset between the charged and neutral velocities 
in Fig. 12c, o? is a numerical artifact. See text.) 
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Fig. 13. — Time evolution of a shock in the WC and SC models. At time t = the shock 
discontinuity is located at 2; = and is propagating in the +z-direction. (a) Velocity profile 
in the WC model a,t t — 6.8 yr. (b) SC model at the same time, (c) WC model at i = 68 yr. 
(d) SC model at the same time. 



-40- 



U = 272 yr 




01234 01234 
z (10^^ cm) z (10^^ cm) 

= 4014 yr 




z (10 cm) z (10 cm) 



Fig. 14. — Same as in Fig. 13, but at later times, (a) WC model at 272 yr. (b) SC model 
at the same time, (c) WC model at 4074 yr. (d) SC model at the same time. 
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Fig. 15. — Same as in Figs. 13 and 14, but at later times, (a) WC model ai t — 5.43 x 10^ 
yr. The shock is still in transition, having the characteristics of a J shock with a magnetic 
precursor, (b) SC model at the same time. By this time, the shock has become a steady 
C shock, (c) WC model at 2.58 x 10^ yr. It is now also a steady C shock by this time. It 
takes much longer to form a C shock in the WC model than in the SC model because the 
coUisional drag forces (which accelerate the neutrals in the precursor) in the WC model are 
less than in the SC model. 



-42 - 




Fig. 16. — Location of the "front" ;i;front(^) of the hydromagnetic signal that propagates 
from the initial shock discontinuity at 2; = and t = for the WC {filled boxes) and SC 
{open boxes) models displayed in Figs. 13 - 15. The location of the front is taken to be the 
point at which the magnetic field strength B first increases perceptibly (~ 0.5%) above its 
preshock (i.e., undisturbed) value. The theoretical prediction of 2;front(0 for the ion-diffusion 
mode (oc t^^^, see eq. [27]) is displayed as the solid hne. The dashed curve corresponds to 
^front = ^disp + V^ms^j which is the prediction for the grain magnctosound mode, displaced by 
the amount ^^disp to coincide with the location of 2;front at the time t ~ ^g~^- 
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Fig. 17. — C shock flow in "intcrmcdiatcly-coupled" models. The initial shock conditions 
and physical parameters are the same as for the WC and SC models shown in Figs. 13 - 
15, except that the IC-1 model has grains with Og = 0.125 //m (Fg = 1.1), while IC-2 has 
Og = 0.09 iim (Fg = 1.3). (a) IC-1 model when it has a steady C shock at t ^ 1.70 x 10^ yr. 
(6) Steady C shock flow in the IC-2 model at t = 1.25 x 10^ yr. 



